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ABSTRACT 


Electric currents induced by non-periodic winds in the 
ionosphere are calculated for the case of incoincidence of 
the Earth's rotational and magnetic axes, based on a real- 
istic model of the upper atmosphere. Main results obtained 
are as follows: (1) Sq-like current systems changing with 

universal time are produced and the intensity of main vor- 
tices is about one tenth of that of the 3q field for a typi- 
cal wind velocity of 10 m/sec. (2) Of the two components of 
non-periodic winds the meridional component is more effective 
in producing current systems. (3) The position of main cur- 
rent vortices is controlled by the distributions of ionos- 
pheric conductivity and of wind velocity and the degree of 
influence of these two is different for different profiles 
of winds. (4) The intensity of main vortices changes with 
longitude in such a manner that the maximum intensity occurs 
in the north and south American zone and the minimum in the 
Asia and Oceania zone, being in agreement with that of the Sq 
field. (5) No remarkable difference can be seen in the results 
obtained by using different coordinate systems. 
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1. INTRODUCTION 

The atmospheric dynamo theory was studied by a number 
of authors (see for example, a brief summary by Maeda . 1966) 
for interpreting the solar and lunar daily geomagnetic varia- 
tions at the early stage, and for discussing dynamical beha- 
vior of the upper atmosphere at the latter. In most of these 
studies only periodic components of ionospheric winds which 
induced the so-called dynamo electric field were considered, 
because it has long been believed that ionospheric winds are 
of only or mainly tidal origin. 

In recent years, however, thermally driven ionospheric 
winds have warmly been discussed, and it has been pointed out 
by Gupta (1967), for example, that the contribution of gravi- 
tational tide to (solar daily) geomagnetic variations would 
be negligible as compared with that of solar radiation. If 
it is so, the winds may have non-periodic components as well 
as periodic ones. In fact the results of observation of iono- 
spheric motions by radio and rockets strongly suggest the 
existence of non-periodic components in ionospheric motions. 
Thus, it may be expected that electric currents might be in- 
duced by such non-periodic winds and may result in geomagnetic 
variations observed on the ground. 
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This problem was considered by Kato (1957) for a zonal 
wind, and by Jones (1963) and DeWitt and Akasofu (1964) for 
a solenoidal wind, and it was concluded by these authors that 
j no total currents could be induced by such a zonal or sole- 

noidal wind. Their discussions were made on some simplifying 
assumptions, especially it was assumed that the electric cur- 
rents flowed in a thin layer in which the wind velocity and 
the electric field did not change with height (a two-dimen- 
sional assumption ) and that a rotational and magnetic axes 
of the Earth were coincident (a coincident-axes assumption ) . 
There still, therefore, remains the possibility that a cur- 
rent system can be produced by non-periodic winds, for the 
case in which one or two of the above assumptions are not 
satisfied to a certain approximation. In fact, the results 
of estimate by Maeda and Matsumoto (1962) and by van Sabben 
(1962) seemed to suggest this possibility. 

Maeda and Matsumoto considered the case in which the 
coincident-axes assumption did not hold, and they obtained 
current systems depending on universal time as well as local 
time, van Sabben discussed the case in which the two-dimen- 
sional assumption was violated, and he showed that certain 
height-dependent zonal or meridional wind distributions would 
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cause Sq-like current systems. These findings seem to be im- 
portant in discussing not only the daily geomagnetic varia- 
tions but also the dynamics of the upper atmosphere. Their 
numerical results, however, seem to be far from any realistic 
one, because the models employed were over-simplified. 

The purpose of this series of papers is to discuss the 
method of mathematical treatment of the problem in a general 
form and then present numerical results based on a realistic 
model of the upper atmosphere. Part I deals with only the 
case of incoincidnet axes. The effect of height-dependent 
wind structure will be considered in Part II where the axes- 
incoincidence is ignored, and the most general case in which 
the above two cases are combined will finally be discussed 
in Part III. 


2. TWO-DIMENSIONAL DYNAMO EQUATIONS 

In the two-dimensional case, we take a hypothetical 

spherical current sheet which has a height-integrated tensor 

conductivity £Kj , then the (height-integrated) current densi- 
— > 

ty J is given by 

? = £K] E t (1) 

where E t is the (height-independent) total electric field. 
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If we consider a quasi- steady state, the effect of self- 
inductance may be ignored, so that the total electric field 

is given by the sum of the electrostatic field E s derived 

— > 

from a potential S and the dynamo electric field E^ (= 

— ► — ► — > 

V x B; V is the wind velocity and B is the geomagnetic field 

intensity) . Since displacement currents can be neglected in 

the present case, the current must be divergent free; i.e., 

— > 

div J = 0. (2) 

Taking spherical coordinates (co-latitude & , longitude A ) 
and local Cartesian coordinates ( x for southward, y for 
eastward) at a point ( 5,A) on the current sheet, we have 
from (1) and (2) 
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and Kj cx , K and K x ^ are components of tensor conductivity 
(see for example, Chapman , 1956). Since oi B y 0 this 
is a partial differential equation of elliptic type. 

The electric conductivities (K^ etc.) are, in general, 
functions of Q , A and t (local time) , but the dynamo field 
does not depend on local time because only non-periodic 
winds are considered here. Thus eq. (3) may be written 


5 0 * d O <='/'•- 


(4) 


where 


A = |3/o<=A(0,A>t) 

B=T/oC=B(0,A» t ) 
c = S /oc = c ( e • X > t) 

D = £ /oC = D ( e > X ' t ) 

and therefore the solution would be the form 

S = S ( 6 , A , t ) . (5) 


Eq. (4) is the basic equation in our problem, to be solved 
under appropriate boundary conditions. 

In actual cases the coefficients. A, B, C and D cannot 
be expressed in any analytical forms but are given by numeri- 
cal values, and therefore the solution is obtained by a method 
of numerical integration. When we carry out calculations. 
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an important question may arise, i.e., which of the geo- 
graphic and geomagnetic coordinates is appropriate for the 
practical treatment ? We cannot briefly answer this ques- 
tion, and this will be discussed in the next section. 

3. COORDINATE SYSTEMS AND BOUNDARY CONDITIONS 
The physical quantities involved in eq. (4) are the 
wind velocity V, the geomagnetic field intensity B, the 
electrical conductivity CKJ as known, and the electrostatic 
potential S as unknown. Of these, the wind velocity would 
be dependent on the geographical (gg) coordinates, and the 
geomagnetic field intensity depends, of course, upon the geo- 
magnetic (gm) coordinates. Thus, if we take the gg coordi- 
— ^ 

nates, B must be transformed; whereas if we take the gm co- 
— ^ 

ordinates, V must be transformed. However, what about the 
electrical conductivity ? The electrical conductivity con- 
tains several quantities; some of them (for instance, the 
electron density or collision frequency) might be dependent 
on the gg coordinates, but the gyrofrequency would be in con- 
formity with the gm coordinates. Thus it seems to be very 
hard to determine any one coordinate system appropriate for 
the conductivity. We know, however, that near the equator 
the conductivity is strongly controlled by the geomagnetic 
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field as well known as the equatorial electrojet phenomenon, 
whereas in high latitudes the gg coordinates might be better 
for the expression of conductivity, because the geomagnetic 
field is nearly vertical there. 

We meet the same difficulty as above for the determi- 
nation of boundary conditions. If we want to solve the 
basic equation in its original form, we must employ one of 
the following two kinds of boundary conditions: One is 

called a Dirichlet problem in which the values of the un- 
known at the boundary are known and the other is called 
a Neuman problem in which the values of the normal derivative 
of the unknown at the boundary are known. Unfortunately, our 
present problem is not so simple. The unknown S is a poten- 
tial, so that we can take an arbitrary value of S at any one 
point. The north pole is usualy taken as this point. However, 
a difficulty again occurs, because the values of all the coe- 
fficients (A, B, C and D) in the equation tend to infinity at 
the pole, so that we cannot set any boundary conditions at 
the pole. Instead, the polar condition might be replaced by 
the values at points close to the pole, and these values may 
be estimated by Taylor expansion of S around the pole. This 
boundary condition would be better satisfied for the gg coor- 
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<3inates rather than the gm coordinates. 

On the other hand, at the equator we may take the 
dition that 

as/ae = 0 

which is based on the assumption that the geomagnetic field 
lines can be regarded as equipotentials, and therefore this 
condition is satisfied for the gm equator not for the gg 
equator. Thus, we again meet a difficulty for the choice 
of coordinates appropriate for boundary conditions. Since 
i® impossible to take any combined coordinate system, cal— 
culations are carried out for each of the gg and gm coordi- 
nates. We would say therefore that the results taking the 
gg coordinates are doubtful near the equator, whereas the 
results taking the gm coordinates are less significant in 
the polar region. 

Finally, for the east-west boundaries we set the condi- 
tion 

S (at A= 0) = s (at A = 27c) 
because of the periodicity expected for S. 

4. COORDINATE TRANSFORMATION 


As discussed above, if we employ any one coordinate sys- 



tem, some quantities must be transformed; i.e.. 


if we take 


the gg coordinates, the geomagnetic field intensity B must 
be transformed; whereas if we take the gm coordinates, the 
wind velocity V must be transformed. The method of trans- 
formation is discussed below. 

4. 1. General . 

A general method for the coordinate transformation of 
any physical quantities as a function of position on a 
sphere may be provided by using the spherical surface 
harmonics (see, Schmidt . 1935). 

Let N and N' in Fig. 1 be the north poles of any two 
coordinate systems (u,t) and (u',t'), we have 


Yl 

P^ (cos u') cos u- 1 1 = Ara/i p n m ( cos cos 

n y yn~o / 

Yl 

(cos u*) sinyttf =2 B m/t P n m ^ cos sin 


mt 


mt 


( 6 ) 


where P m are the semi-normalized associated Legendre 
n 

functions introduced bv Schmidt (1935), and A n m and B n m are 
given by 


V" V» " (a ™ v> <A + B)/2 

■y ■ <a ™ v )>5 (A ' B)/2 

A .a-cr/-^0-cT£!^>] 
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B - (1 + cf /- [ (1 - cl* amp n( c > 1 

L d c m J 

c = cos e 
s = sin e 

a nm = £m <n ~ m )!/( n + ">) ! 

£.“ 1. £,-€»“ - 2 

and P n (c) = Le 9 en< 3 re functions. 

These formulae are the most general expression for the 
transformation of spherical coordinates, and they are very 
useful for the present problem. 

4. 2. Case where the gg coordinates are taken . 

If we take 

u = 6 9 u ' = (S) 

t = A-^O; t ' = 7t — A 

e = B 0 , 

where ( 9 , A ) are the gg coordinates of point P (see Fig. 1) . 
((E), -A.) are its gm coordinates, and ( 9 0 ,\ 0 ) are the gg coor- 
dinates of the geomagnetic pole N' , then we have from (6) 

yi 

(cos © ) coSy«( 7 r-A ) = 2 < cos 0 ) cos m ( A - Ao) 

> "' 0 (7) 

It 

P^(cos©) sinyjt(7C~A) =2 B m ^ p /U cos 0 ) Sin m(A-Ao) 

w = o “/ ' 

The geomagnetic field intensity can, in general, be trans- 
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formed into the gg coordinates by these formulae. 

As an example, if we approximate the geomagnetic field 

— > 

by a dipole field of magnetic moment M, we have 

B z = -(2M/a 3 ) P 1 °(cos@) 

for the vertical intensity. Application of (7) to P^(cos©) 
gives 

P 1 °(cos^)) «P 1 °(c) P 1 °(cos^) 

+ P-^Cc) P^ (cos 0 ) cos(A-Ao) 

= cos $0 cos# + sin^ sin 0 cos (A “Ac) 

it follows that 

B z =-(2M/a 3 ) £cos0 o cos# + sin# o sin# cos ( A - 

This is the well known expression for the vertical intensity 
of the goemagnetic field transformed into the gg coordinates. 

4. 3. Case where the gm coordinates are taken . 

The inverse transformation to the above can be obtained 
by changing 

t t ' , u ( — t u ' 

and again by taking 

u = 6 j u ' = ©, 

t = A - A 0 , t 1 = 'K - A- 

e " Qo • 
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The results are 

Tt 

p n ( c °s 6 ) cos (A-Ao) A^P m (cos<g)) cos mt'TT-.A) 

VK«0 7 11 

n (8) 

P n (cos G) sin =2^B n ^P n m (cos(g)) sinm(TT-A) 

These formulae can be used for the transformation of the wind 
velocity into the gm coordinates in the general case. 

If we take = 0 (for non-periodic winds) we have 

P n ( cos ^ ) = ^n" 1 ^ p n m ( cos @) cos m(TC-A) 

-Wt' ° 

Taking n = 1 for example, 

P^tcosfl) = P 1 < ^(c) P 1 °(cos@) 

+ PjV) P 1 1 (cos@) cos('ft'-vl) 

i .e . , 

cos 0 = cos0 o cos(ff} - sin# o sin(S) cos j\ 

and using the relation 

sin 9 - (1 - cos^ 

we have the expression of 2 sin$ cos^ in the gm coordinates. 
This form will be used as an example for the longitudinal 
distribution of the meridional component of non-periodic winds. 

5. CONDUCTIVITY AND WIND PROFILES 
In order to solve the basic equation (4) , we need to 
have the numerical values of the coefficients A, B, C and D, 
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where A, B and C are functions of conductivity only, whereas 

— > 

D includes the wind velocity V. 

The distribution and variation of conductivity are ob- 
tained as follows: Atmospheric parameters, such as tempera- 

ture, density, molecular weight and so on, are taken from 
the U. S. Standard Atmosphere (1962). The distribution of 
electron density with height and latitude at noon are esti- 
mated, for a period of high sunspot activity, from the re- 
sults of recent radio and rocket observations. By using 
these atmospheric and ionospheric models, the height-integ- 
rated conductivities K xx , and are calculated and 

the results are shown in Fig, 2 . Since we do not yet have 
any detailed information about the time-dependent structure 
of the lower ionosphere, especially about the height distri- 
bution of electron density, the following formula for the 
variation of conductivity (Chapman and Bartels . 1940) is used: 

3 s 

K = K 0 (1 + Z k s cos > 
s- I 

where 9C is the zenith distance of the sun. The values of 
the coefficients k g are estimated as follows, by taking 
into account the latitudinal distribution of noon conductivity 
and the height distribution of electron density observed by 
rockets at night (e.g., Bourdeau, 1963): 
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k x = 3.0, k 2 = 3.0, k 3 = 0.9 

This gives the night-time conductivity of about one fifties 
of the daytime one at middle and low latitudes. 

Although we have now so many results of observations of 
winds at heights lower than 100 km (see for example, a review 
by Murgatrovd . 1957) , very few are known about the wind struc- 
ture in regions higher than 100 km. This seems to be due to 
the reason that the method of wind measurements taken at the 
lower atmosphere cannot be applicable to the upper atmosphere 
because of low density, and that the wind structure at heights 
of interest (say, around 120 km) is very complicated perhaps 
because of the interaction of the neutral and ionized compo- 
nents of atmospheric gases. Thus, unfortunately, we cannot 
as yet set any reliable profiles of non-periodic winds. We 
shall employ here, therefore, only hypothetical wind profiles 
as shown in Table 1 and illustrated in Fig. 3 . which might be 
deduced from informations about ionospheric drifts and also 
wind structures in the lower atmosphere. Cases (AZ) and (AM) 
are chosen to be similar to those adopted by van Sabben (1962) 
for the convenience of comparison. 
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Table 1. Supposed profiles of non-periodic winds. 



Meridional (Southward positive) 

Case (AM) 

v x = - v 0 sin ( 2 9 ) 

Case (BM) 

V x = - V Q sin (40) 


Zonal (Eastward positive) 

Case (AZ) 

V y = - V Q (2.7 sin# - 4.0 sin 3 0 ) 

Case (BZ) 

V y = - V 0 sin (40) 


Note: Q is co- latitude. 

V is taken to be 10 m/sec in this paper. 


6. CALCULATION OF THE DYNAMO FIELD 
In the case where the gg coordinates are taken, the velo- 
city distribution can be used in its original form, and we 
have the following expressions for the dynamo fields: 


where 


E dx = -C F V y 


E , = C F V 

dy x 


, 3 

C = 2 M/a , 

F = cos Q 0 cos 9 + sin sin Q cos ( A 


which has been obtained in section 4.2. 
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In the case where the gm coordinates are taken, the 
velocity distribution must be transformed. This transform- 
ation can be made by using the method presented in section 
4.3., and we have the following results: 


where 


E dx 

- V o C 

F, F 
£ 1 r 5 

E , 

= - V 

C F. 

dy 

o 

3 

E, 

= 4 V 

C F, 

dx 

o 

1 

E, 

dy 

II 

1 

o <; 

o 

*1 




2 V C F. 
o 2 

2 V C F, 
o 1 

4 V C P n 
o 2 

4 V CP, 
o 1 



P 3 P 4 \ 
P 2 P 4 J 


for Case (AZ) 


for Case (BZ) 


for Case (AM) 


for Case (BM) 


F^ = cos(§)(cos0 o sin(£) + sin Q o cos $ cosy).) 

F_ = cos£ cos @ - sin# sin @ cos A 
F 3 = sin Q o cos <g) sin A 
F 4 = 2 (P 2 ) 2 - 1 
p 5 = 4 ( p 2 ) 2 “ 1.3 

7. METHOD OF NUMERICAL SOLUTIONS 
There are two methods for numerical solutions of eq. (4) . 


One deals with the equation in its original form, and the 
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other is to modify the original equation to ordinary dif- 
ferential equations by expanding the potential S (unknown) 
into Fourier series. It seems that the latter is, in general, 
more easy to obtain solutions than the former. In this case, 
however, we must interrupt the series at the first several 
terms for practical calculations, so that significant error 
might be introduced. For this reason, we adopt here the 
former method, that is, eq. (4) has been modified to a finite- 
difference form and solved by Liebmann's iteration method on 
a IBM system/360 computer. We do not want to describe details 
for the method employed, but give only essential points below. 

Returning to eq. (4) , if we use the relation 

t = ut + A* 

among the local time (t), the universal time (UT) and the lon- 
gitude (A), the coefficients A, B, C and D are regarded as 
functions of space coordinates ( $ , A. ) only, at a particular 
UT. Thus, the numerical integration can be made only for space 
coordinates, and we have a set of solutions for different UT's. 
The mesh has been taken to be 2.5° for Q , and 7.5° for A. # 
and solutions have, in most cases, been converged within limits 
of error for several minutes. This would be due to the reason 
that our basic equation is, fortunately, of elliptic type. 
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Concerning the polar boundary condition, we have tried 
the case in which all the boundary values are taken to be 
zero, and found that no significant difference is seen in 
solution and also in time of convergence, as compared with 
the case in which the boundary values estimated by Taylor 
expansion are used. This is due to the fact that the boun- 
dary values close to the north pole are very small, because 
we have assumed that S = 0 at the pole. 

8. RESULTS AND DISCUSSIONS 

The distributions of the electrostatic potential (S) 

ll 

and the current functions (R) at UT = 0 calculated for dif- 
ferent wind profiles, as is shown in Table 1 , are illustrated 
in Figs. 4A to 4D . These results are obtained by using the 
geomagnetic coordinates and by taking V q = 10 m/sec. It is 
found from these figures that: 1) For all of these four cases, 
current vortices are produced and their patterns are similar 
to those of Sq, namely, they have counter-clockwise flow of 
currents and much enhanced intensity in daytime. 2) The in- 
tensities of daytime current cortices are 

6,000 amperes for Case (AZ) 

16,000 amperes for Case (AM) 
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4.000 amperes for Case (BZ) 

8.000 amperes for Case (BM) 

so that the meridional wind is more effective in producing 
current systems than the zonal one, and these intensities are 
about one tenth of those of the Sq vortices. 3) In most cases 
clockwise and weak current vortices are produced in high lati- 
tudes. 


The cases whe r e the <iuual and »>«. iuiuua 1 wiuuo aic com- 
bined are also calculated and we have Fig. 5 for R. The in- 
tensity of main current vortices is 20,000 amperes for Case 
(A) and 10,000 amperes for Case (B) . It is seen that the 
center of current vortices is much closer to the equator for 
Case (A) than for Case (B) , and it may therefore be said that 
the distribution of conductivity had a leading effect for 
Case (A) in determining the position of current vortices, where- 
as the distribution of wind velocity played more important role 
for Case (B) . In this connection, the Sq like current systems 
obtained by van Sabben (1962) might be much distorted, if the 
latitudinal and local time variations of conductivity are taken 
into account, even in daylight hours. 


In order to see the effect of different coordinate systems 
on the results, similar calculations are made by using the geo- 
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graphic coordinates for Case (AM) and Case (BZ), as an exam- 
ple, and we have the results as shown in Fig. 6 , Comparison 
of Figs. 4 and 6 for corresponding cases shows that no re- 
markable difference can be seen in the intensity and also in 
the general pattern of current systems, though there are 
slight differences in the shape of current vortices. 

All the figures presented above are the results obtained 
for UT = 0* 1 . An example for different UT's is shown in Fig. 
2 , where © shows the position of noon. It is found that 
both the shape and the intensity of current systems are chan- 
ging with UT, and this may correspond to the UT variation in 
the Sq field. In order to see this point more clearly, we 
separate these changing current systems in such a way as 
follows: 

First, if we fix the local time and average over UT's, 
then a current system as shown in Fig. 8A is obtained. Since 
this system depends only on local time, it may be called the 
local time part and this part would have a main contribution 
to the Sq current system. 

Next, if we fix the longitude and average over UT's and 
subtract the non-periodic component of R, because this part 
has already been included in the local time part, we have the 
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results as shown in Fig. 8B . This system depends only on 
longitude and therefore it may be called the longitudinal 
part . Since this part gives a longitudinal inequality in 
current systems, it is expected that the maximum intensity 
of current vortices, and therefore the madimum magnetic effect, 
occurs around j\. (geomagnetic longitude) = 0°, i.e., the north 
and south American zone, and the minimum around A.- 180°, 
i.e., the Asia and Oceania zone. This tendency is in good 
agreement with that as seen in the Sq field (see for example. 
Price and Wilkins . 1963; Matsushita and Maeda . 1965) and also 
in the equatorial electrojet (Sugiura and Cain . 1966). It may 
therefore be said that the longitudinal inequality in the Sq 
field is caused partly by the incoincidence of the Earth's 
rotational and magnetic axes, and partly by the distribution 
of ionospheric conductivity due perhaps to the longitudinal 
inequality in the geomagnetic field intensity. 

Finally, if we take the non-periodic part at each UT and 

subtract the average over UT's, we have the universal time 

part depending only on universal time as is shown in Fig. 8C . 

It is expected from this figure that the current vortices would 

ll 

be most enhanced around UT - 20 and least enhanced around UT 
= 9 . Note, however, that these three parts cannot be indepen- 
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dent because of the relation, t = UT + _A , where t is geomag- 
netic local time and UT is geomagnetic universal time defined 
by refering to the meridian = 0°, because the geomagnetic 
coordinates are taken. 

9. CONCLUSIONS 

From the foregoing detailed calculations of electric 
currents induced by non-periodic winds in the ionosphere, 
when the incoincidence of the rotational and magnetic axes 
of the Earth is taken into account, the following conclusions 
may be drawn: 

1) Sq-like current systems can be produced by non- 
periodic components of ionospheric winds as an effect of the 
axes incoincidence, even for the two dimensional case of the 
dynamo theory. 

2) When a typical wind velocity of 10 m/sec is taken, the 
induced current intensity is about ten thousand volts, being 
about one tenth of that of the Sq field. 

3) Of the two components of non-periodic winds, the mer- 
idional component is more (two or three times) effective in 
producing current systems. 

4) The position of main current vortices is controlled 
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by the distributions of ionospheric conductivity and of wind 
velocity, and the degree of incluence of these two is diffe- 
rent for different wind profiles. 

5) No remarkable difference can be seen in the results 
obtained by using different (geographic and geomagnetic) coor- 
dinate systems. 

6) The intensity and the pattern of current systems change 
with longitude and universal time. The maximum intensity 
occurs around 0° in geomagnetic longitude (the north and south 

American zone) and 20^ in geomagnetic universal time, and the 

o h 

minimum around 180 (the Asia and Oceania zone) and 9 , res- 
pectively. 

Since the intensity of current vortices is proportional 
to the velocity of winds, it would be said that the effect of 
axes incoincidence cannot be negligible when much stronger 
winds exist. 
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FIGURE CAPTIONS 


: 


Fig. 1. Showing relation of any two spherical coordinate 
systems, (u,t) and (u',t'). 

Fig. 2. Latitudinal distribution of height-integrated noon 
conductivities, K^, Kyy and K X y* 

Fig. 3. Illustration of the four profiles of non-periodic 
winds as is shown in Table 1. 

Fig. 4A. Distribution of the electrostatic potential (S) 

and the current function (R) for Case (AZ) , where 
electric currents flow in the direction indicated 
by arrows along the equal R (strem) lines. 

Fig. 4B. Distribution of S and R for Case (AM). 

Fig. 4C. Distribution of S and R for Case (BZ) . 

Fig. 4D. Distribution of S and R for Case (BM) . 

Fig. 5. Distribution of the current function (R) for com- 

bined winds; Case (A) is obtained by a combination 
of profiles (AZ) and (AM) , and Case (B) by a com- 
bination of profiles (BZ) and (BM) . 

Fig. 6. Distribution of the current function (R) calculated 
by using the geographic coordinates for wind pro- 
files (AM) and (BZ) . 

Fig. 7. Variation of current systems with universal time, 
where (•) shows the position of noon. 

Fig. 8A. Local time part of the current systems as shown in 
Fig. 7. 

Fig. 8B. Longitudinal part of the current systems as shown 
in Fig. 7. 

Fig. 8C. Universal time part of the current systems as shown 
in Fig. 7. 
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